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Xt = e" 5t ^ Xo + J e^-du 

where £ s is a Levy process. Our primal goal is to estimate the characteristics 
of the Levy process £ from the low-frequency observations of the process 
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1. Introduction 


Let (£t)t>o be a Levy process with a Levy triplet (p,, The main object 


of our study is the so-called generalized Ornstein-Uhlenbeck (GOU) process 
defined as 



(1) 


The GOU processes have recently got much attention in the literature. A com¬ 
prehensive study of the GOU processes and an extended list of references can 
be found in the theses of Behrne [2], where, in particular, it is shown that X t 
satisfies the following SDE: 



0<s<t 


The popularity of GOU processes is related to the fact they appear to be useful 
in several applications. For instance, the process (1) determines the volatility 
process in the COGARCH (COntinious Generalized AutoRegressive Condition¬ 
ally Heteroscedastic) model introduced in Kliippelberg et al. [18]. One impor¬ 
tant result from the theory of GOU processes is that, under some conditions, 
the process (1) is stationary with invariant stationary distribution given by the 
distribution of the following exponential functional of £ : 



( 2 ) 


In fact, the properties of the functional (2) have been widely studied in the 
literature and we refer to the survey by Bertoin and Yor [7] for a theoretical 
background of the exponential functionals. In particular, it is known that the 
Mellin transform of the density tt of exponential functional, 



satisfies the following recursive formula 


M(z) = ^ M(z + 1), 
z 


( 3 ) 


where </>(z) is a Laplace exponent of the process £, i.e., cj)(z) := — logE [e Z?1 
and complex z is taken from the strip 


T := jz G R : 0 < Re(z) < 0 j with 0 := sup {x > 0 : Efe"^ 1 ] < l} . (4) 
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The recursive formula (3) first appeared for real z in the paper by Maulik and 
Zwart [25]. The validity of (3) for complex z was recently shown by Kuznetsov, 
Pardo and Savov [20]. If is a subordinator, the parameter 9 is equal to infinity. 
Let us note that the functional appeared in such application areas as finance 
(see, e.g. the monograph by Yor [33]), carousel systems (see Litvak and Adan 
[23], Litvak and van Zwet [24]), self-similar fragmentations (see Bertoin and Yor 
[7]), and information transmission problems (especially TCP/IP protocol, see 
Guillemin, Robert and Zwart [13]). For the detailed discussion on the physical 
interpretations, we refer to Comtet, Monthus and Yor [10] and the dissertation 
by Monthus [26]. 

In this paper, we mainly focus on the case when £ is a subordinator with finite 
Levy measure. In terms of the Levy triplet (/z,cr 2 ,^), this means that /x > 0, 
a = 0, v(M-) = 0 and moreover A := v{]R + ) < oo. Suppose that the process 
(1) is observed at equidistant time points 0 = to < t\ < ... < t n . Since under 
some mild assumptions the process is stationary and the invariant distribution is 
given by the distribution of the exponential functional A (see Fasen, [12]), we 
assume that X to ,..., X tri are also distributed as A. Our main goal is statistical 
inference on the Levy triplet (/x, er 2 , v) based on the observations X to , ..., X tn . 
More precisely, we will pursue the following two aims: (1) estimation of the drift 
term /x and the intensity parameter A; (2) estimation of the Levy measure v. 

To the best of our knowledge, the statistical inference for GOU processes of 
the form (1) from their low-frequency observations has not been yet studied in 
the literature. In fact the resulting statistical problem is quite challenging and 
needs careful treatment. Indeed, the only connection between the stationary 
distribution of a GOU process, which can be estimated from the data, and the 
parameters of the underlying Levy process is given by the recurrent relation 
(3) which is rather implicit. The main idea of our procedure for estimating the 
parameters of the process £ can be described as follows. First, by making use 
of (3), we estimate the Laplace exponent cf>(z) at the points z = u° + iv G T, 
where u° > 0 is fixed and v varies on the equidistant grid between eV n and V n 
(with s > 0 and V n —► oo as n —> oo) Afterwards, we use the representation 

(f>(u° +\v) = A + n (u° + it;) — T[v\(—v), v € ]R, (5) 

where v(dx) := e~ u X u(dx), and T[v\(v) stands for the Fourier transform of the 
measure v , i.e., F[v\(v) := j R e lvx u(dx). Since F[v\(v) —> 0 as v —> oo by the 
Riemann-Lebesgue lemma, upon taking real and imaginary parts of the left and 
right hand sides of (5), we are able to consequently estimate the parameters 
/x and A. With no doubt, the second aim, a complete recovering of the Levy 
measure v 1 is the most difficult task. Since the estimates of the parameters [i 
and A are already obtained, we can estimate by (5) the Fourier transform J-[v\[y) 
of v for v from [— V n , V n ]. The last step of this procedure, the estimation of the 
Levy measure v, is based on the regularised inverse Fourier transform formula. 

The above estimation algorithm bears some similarity to the spectral estima¬ 
tion algorithm introduced by Belomestny and Reiss [4] , [5]. Let us also mention 
that the problem of statistical inference for Levy processes (or some their gener¬ 
alizations) observed at low frequency was the subject of many studies, see, e.g. 
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Neumann and Reifi [27], Reifi [28], Kappus [16], Trabs [31] and Jongbloed et al. 
[15]. Note that the last reference deals with the Levy-driven Ornstein-Uhlenbeck 
processes, which are not of the form (1). 

The paper is organized as follows. In the next section, we formulate our main 
assumptions and give some examples. In Section 3, the main estimation algo¬ 
rithm is presented and discussed in details. Next, we analyze the convergence 
rates of the proposed algorithms in Section 4 and provide some numerical exam¬ 
ples in Section 5. The proofs of our theoretical results are collected in Section 6. 

2. Main setup 


In this article, we study the class of subordinators with finite Levy measures as 
possible choice for the Levy process (£ t ). In terms of the Levy triplet (/x, A, zz), 
this means that 


f t 1 > 0, <7 = 0, 

| v(]R-) =0, A := zz(IR + ) < oo. 


(6) 

A detailed discussion of the subordination theory as well as various examples of 
such processes (Gamma, Poisson, tempered stable, inverse Gaussian, Meixner 
processes, etc.), are given in [1], [6], [11], [29], [30]. Note that in the case of 
subordinators, the truncation function in the Levy-Khinchine formula can be 
omitted, that is, the characteristic exponent of £ is equal to 

/»00 

ip(z) = logE [e 12 ^ 1 ] = i [iz + / ( e lzx — l) v(dx). 

Jo 


(7) 

Later on, we also need the Laplace exponent of £, which is defined 

as 


(f(z) := — logE [e -2 ^ 1 ] = —if (iz). 



Under the assumption (6) the Laplace exponent <p(-) is given by 



r- OO 

cfi(z) = tiz+ (1 — e~ zu ) v{du) 

Jo 


(8) 

/»oo 

= nz + z / e~ zu v ((u, +oo)) du. 

Jo 


(9) 


Let us summarise the main properties of the functional A = J 0 °° e dt in this 
case. 


Proposition 2.1. The random variable A ^ admits a bounded density 7r and 
fulfills E [Aj^ 1 ] < oo for all s > 0. If n > 0, then 0 < Aoa < l//i a.s. Moreover, 
the following relation holds for Re[z] > 0, 

m = 112+ 1 (1 “ e "“ w<fc) ’ 

/»oo 

= A +/j,z — I e~ zx v{dx), (10) 

J o 

where A i(z) is the Mellin transform ofn. 
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3. Estimation of the Levy triplet 

In the sequel, we suppose that we are given by the observations X to , X tl ,..., X trl 
of the process (1) at the equidistant time points 0 = t 0 < t± < ... < t n with 
tj = j ■ A for some A > 0. Assuming that the process X t is stationary (see [2], 
[12]), we get that the random variables X^ := X tk , k = 1 ,n, have all the 
same distribution, which coincides with the distribution of A^. 

The first step of our estimation procedure consists in the estimation the 
Laplace exponent (f>(z) for z = u° + iu, where u° > 0 is fixed and v varies. An 
estimator of 4>(z) can be obtained from the recursive formula (3) for the Mellin 
transform of n. In [9], this formula is proved for real positive z such that cf>(z ) > 0 
and M(z + 1) < oo. The case of complex 2 is considered in [20], where one can 
also find some generalizations of the formula (3) to the integrals with respect 
to Brownian motion with drift. In particular, applying Theorem 2 from [20], we 
get that (3) holds for any z 6 T. In the situation when (£ t ) is a subordinator, 
the set T coincides with the positive half-plane (equivalently, the parameter 6 
is equal to infinity) due to 

E [e -2 ^ 1 ] = —4>(x) = —fix — x f e~ xu v ((u, +oo)) du < 0, V x > 0. 

Jm.+ 

Motivated by (10), we first estimate the Mellin transform A i(z) via its em¬ 
pirical counterpart 


1 

M n (z) := - VX 


Z- 1 

k 


(ii) 


and then define an estimate of the Laplace exponent <j>(z) by 


Y n (z) 


M n (z) 
A4 n (z + 1) 


( 12 ) 


If the sequence X \,..., X n has some mixing properties, then we can expect that 
Y n (z) —> 4>{z) in probability. 


3.1. Estimation of A and /x 

Under our assumptions, the Laplace exponent of the Levy process (£t) can be 
represented in the form: 

(f>(u° + iv) = A + /x • (u° + iv) — J r [i>](—v), v £ M, (13) 

where A := f M v{dx) and v(dx) := e~ u ° x v(dx). The general idea of the pro¬ 
cedure described below is to estimate the Laplace exponent <?!>(•) at the points 
z = u° + iv with v £ R and then use the relation (13) for the estimation of 
parameters. By the Riemann-Lebesque lemma, F[v}{—v) —> 0 as v —> +oo (see, 
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e.g., [17]) and we conclude from (13) that cj>(u° + iv) is approximately (at least 
for large v) a linear function in v with the slope /i and the intercept term A. 
This observation suggests that a properly weighted least-squares approach can 
be applied to estimate /* and A. Let V n be a sequence of positive real num¬ 
bers and w(-) be a nonnegative weight function supported on [0,1]. Define a 
scaled weight function w n (v) = V~^w(v/V n ) and introduce the estimators of 
the parameters A and /i as the solution of the following optimization problem: 


(^n i /*n ) ■ 


f°° 2 

argmin / w n {v) \Y n (u° + iv) — ^ ■ (u° + iv) — A|“ dv 
(A,/i) J o 


arg mm 
(a,m) * 


w, 


t (u) {| Im [y n (u° + iv)] - fj,vf 


+ |Re [Y n (u° + iu)] — A — /j,m°| 2 | dv 


with Y n (z) defined in (12). The above optimisation problem admits an explicit 
solution given by 


Mn — 


— 


/ 0 °° u;„(u)Im [Y n (u° + in)] dv 
/ 0 °° vw n (v) dv 

f 0 X w n (v )Re [Y n (u° + iu)] dv 


Wn,n{v )Im [Y n (u° + iv)] dv 


f 0 x w n (v)dv 


II n U 


/»oo 

/ w\, n (v)Re [Y n (u° + iv)] dv - n n u° 
Jo 


with w^ n {v) := V n 2 w fl (v/V n ) and try„(v) := V n 1 w x (v/V n ), where 


.(•) = 


-i 


(•), wx(-) = C 0 w w(-), c ifW = / v l w(v) dv, * = 0,1. 


Taking into account the dehnition of the weight function w n {•), we get also some 
equivalent representations of the estimators /i„ and \ n 


l^n 


A 


n 


arg mm 

A* 



u + iaV n 


)] 



2 

da 


arg mm 
A 



- « - 


2 

da. 


In practice, we need to replace the above integrals by sums. To this end, let 
the numbers ay,..., aM constitute an equidistant grid on the set [e, 1] for some 
e > 0. We estimate the Mellin transform A4(z) for all z £ {u° + ia rn V n , m = 
1,..., M} and z £ {u° — 1 + ia m V n , rn = 1,..., M} and so get the estimates of 
the Laplace exponent at the discrete points z = u° + ia m V n (see above). Now 
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we define an estimate of the parameter /i via 




arg min E ^{pm) (p H"~ fJ'&mVn 

^ m— 1 

1 ^^pnri)Ot-m Im[E 1 (ll -\- ic^ m V^,)] 

V n ■ Em=l M, ( a m)“i 


(14) 

(15) 


Afterwards, we estimate the parameter A by 


A 


n 


M 

arg min E w(a m ) (Re[y„ {lL + io^m^n)] ^ 


m=l 


Em=i w{a m )Re[Y n (u° + ia m 14 )] 

EEi^(«-) 


Hull . 


(16) 

(17) 


The whole algorithm is described below. 
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Algorithm 1: Estimation of A and /i 

Data: n observations Xi,..., X n of the GOU process ( X t ) observed 
at equidistant grid j ■ A, j = 1,..., n. 

Initiate: Fix V n —> oo, e £ (0,1) and u° > —1. 

Set atj = e + j ■ (1 — e) /M, j = 1,..., M. 

Fix a function w(-) > 0 supported on [c, 1]. 

Denote v m ,n ■= a m V n . 

Algorithm: 

1. Estimate the Mellin transform A i(z) := E [A^ 1 ] 
for z £ {u° + Wm^n, 1 + u° + iu m ,„, to = 1,..., M} via 


Mn(z) = -J2* 


Z- 1 

k 


k =1 


2. Estimate the Laplace exponent <fi(z) := — logE [e 2 ^ 1 ] 
at the points z £ {u° + iu mjn , to = 1,..., M} by 


Y n (z) = z 


M n (z) 
■M n (z + 1) 


3. Estimate y by 


Em =1 ™{oim)Otm lm\Y n (u° + 


Mn • — 


Vn ■ 


4. Estimate A by 



M 

m— 1 






3.S. Estimation of the Levy measure v 

As a result of Algorithm 1, we obtain the estimates y n and A„ of the parameters 
/i and A, respectively. Based on (13), we first define an estimate for the Fourier 
transform of v via 

P[v\(-v) = -Y n (u° + m) + y n ■ (u° + iv) + X n . (18) 

Next we estimate the measure v by a regularised Fourier inversion formula 

Vn(x) = — [ e lvx f[v](-v)JC(-v/Vn)dv, (19) 

^ Jm 
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where /C is a regularizing symmetric kernel supported on [—1,1], Note that with 
a slight abuse of notation, we use v also for the density of the Levy measure, 
and v n for an estimate of this density. In what follows, we also use the notation 
v n = e~ u x v n . The formal description of the algorithm is given below. 

Algorithm 2: Estimation of v 

Data: n observations Xi ,..., X n of the GOU process (Xf) observed 
at equidistant grid points j ■ A, j = 1,..., n. 

Initiate: Fix V n —> oo and u° > — 1. 

Set a m = — 1 + 2 ■ j/M, m = 0,..., M. 

Fix a regularizing kernel K, supported on [—1,1]. 

Denote . oi rn \^ n . 

Algorithm: 

1-2 The first two steps coincide with ones of Algorithm 1. 

3. Estimate J r [P](— v m ^ n ) for D(dx) = e~ u X v{dx) by 

^ m,n ) — X n (u -f- Wm^n) fl n ' (U -{- + A n 

for m = 0,..., M. 

4. Estimate v by 

1 M 

K{X) = ■ (1 + M) E 

1 ' 771—0 


Remark 3.1. It is a worth mentioning that the estimation Algorithms 1 and 2 
can be applied to a more general situation when 

= /it+ r t , (20) 

where the process r t is a difference between two subordinators, i.e., r t = r t + +rf, 
and t + and t~ are the processes of finite variation with Levy measures v + and 
v~ concentrated on 1R+ and 1R _, respectively. In fact, in this case, the formula 
(13) still holds with 

v(dx) = E{x > 0 }v + (dx) + E{x < 0 }v~{dx). 

Therefore, the consequent estimation of p, A and the Fourier transform of the 
measure e~ u X v( y dx), as well as the estimation of v are still possible. 

4. Convergence 

In order to analyse the convergence properties of the estimates p n , X n and v n 
we need to further specify the class of Levy processes (£t). 
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Deflniton 4.1. For s £ NU {0} and R > 0, let Q(s,R) denote the set of all 
Levy triplets (p, 0, v), such that v is supported on R + and 

max|^(R + ),y |u| 2s |J r [f'](t’)| 2 < R, (21) 


where v(dx) = e u ° x v(dx). 

Note that if (21) holds, then v is s-times (weakly) differentiable with 


IH S) L<^ jjv\ s \R[i?\(-v)\dv<oc. (22) 

It turns out that the convergence rates of the estimates p n , \ n and v n crucially 
depend on the asymptotic behaviour of the Mellin transform of In order 
to specify this behaviour, let us fix some u° > 0 and introduce two classes of 
probability densities: 


V{p,L) 


p : lirninf [|z;|^ |A 4\p](u° + h>)|] > L 

|d|—>- oo 


(23) 


£{a, L) 


< p : lim inf 

| l>|—>-00 


e“l”l \M[p](u° + iu)| 



(24) 


where a, /3 G 1R, L > 0 and for any density p, A 4\p\ stands for the Mellin 
transform of p. Before we formulate the main convergence results, let us look at 
some examples. 

Example 4.2. Consider the class of Levy processes with p, = 0, a = 0 and the 
Levy density v of the form 


N 


’) = £ 


3 = 1 Lfc=l 


£ 


Qjk'E 


k -1 


e~ PiX ■ I{x > 0} 


with N, rrij G N, pj > 0, g^ > 0. First note that the assumption (6) obviously 
holds. Let us now check (21). We can apply the well-known Erdelyi lemma to 
derive 

[ x k ~ 1 f(x)e lvx dx x civ -fc , v —> oo 


• 1R-i 


for any exponentially decaying and smooth function / on 1R + , and some complex 
Ci depending on /. Therefore, we conclude that 




N rrij 


/ xk 1 fj(x)e lvx dx 

4 — 1 l - 1 J IR.+ 


C 2 V 


-k* 


3=1 k =1 

where fj(x) = e~ ( ' Pj+u ' >x , k* := argmin{3 j : ajk ^ 0} , 

k 
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where > 0 depends on u°. Hence for any s < k* — 1, the condition (21) holds 
for some R > 0. Furthermore, taking into account the asymptotic behaviour of 
the Gamma function (see, e.g., formula 8.328 from [14]): 


|T(u + m)| = exp j — — v + (w — lnu j • V2n (1 + o(l)), 


(25) 


we derive 


N 


K 


| M[u° + if) | x VZnA 1 u exp < —-v + «° ~ x + ^ + E R- e (C?) 


i=i 


3 =i 


where (i,... ,(k are the roots of the equation 

N rn 


g jk (k- 1 )! =a _ 

^ik=i (PJ+*)* 


EE : 


see [19]. Therefore, for any u° > 1/2, we conclude that 7r £ £(ir/2,L) with any 
L > 0. 

Example 4.3. Next, we provide an example of a Levy process with = 
f 0 dt having a density from V((3, L). Consider a subordinator T with drift 
fi > 0 and the Levy density 


v(x) = a6exp{— bx} I{x > 0}, a, b > 0. 


The exponential functional Aqc of the process (£ t ) has a density of the form 
tt(x) = CiX b (l — /{0 < x < 1 //x} 


with some C\ > 0, see [9]. In other words, A <*, has the same distribution as 
£/p, where the r.v. £ has the Beta distribution with parameters a = b + 1 and 
/3 = a/fi = A//u. The Mellin transform of the function n(x) in the half-plane 
Re(s) > —a is hence given by 


M(z) 


e [e- 1 ] 


1 B(z + a - l,/3) 

M* -1 B(a, (3) 

T{a + /3) 1 r(z + a-l) 

r(a) /U _1 T(z + a + /3 — 1 ) 


Using (25), we conclude that the Mellin transform of Aao has a polynomial decay 
in this case. More precisely, 


\M(u°+w)\^L-\v\-W with L = M -“°+ i r(A ^+y 1} 


as |u| —>■ oo and therefore n £ V(\ //z,L). 
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Let us now formulate the main result concerning the convergence of the es¬ 
timates fin and A n . 

Theorem 4.4 (upper bounds for p n and A n ). Let (£ t ) be a Levy process 
with a triplet from Q(s,R). Suppose that the sequence Xq, X^,..., X n is a- 
mixing and strictly stationary. Denote the a-mixing coefficients of the sequence 
X 0 ,X U ...,X n by a(s). 

(i) Assume that the density it of belongs to V(P, L) with some (3 £ JR. and 
L > 0, and moreover 


a(j) < e , j £ N, for some a* > 0. 


(26) 


Then the quadratic risks of the estimates /j, n and A„, under the choice 


V n = n 1 /( 2 / 3 + 2s + 3 ) j satisfy the following asymptotic relations 


E \\p n — p\ 2 \ < n - 2 (s+ 2 )/( 2 / 3 + 2s +3) log(n) 


and 


E |A„-A| 2 < jt, _ 2 ( s + 1 )/( 2 ^+ 2s + 3 ) log(n), 


as n —> oo. 


(ii) If n £ £(a,L) and 


a (J) ~ J a > JSN, for some a* > 2, 


(27) 


then the choice 



leads to the rates 


E \p n -p\ 2 <log 2(s+2) («), 
E [|A n — A| 2 j <log- 2(s+1) (n). 


Proof. Proof is given in Section 6.1. 


□ 


In a similar way, we can establish the upper bounds for the risk of D n . In the 
theorem formulated below, the quality of the estimate D n is measured in terms 
of the mean integrated squared error (MISE) 



Theorem 4.5 (upper bounds for u n ). Let the assumptions of Theorem 4-f be 
fulfilled and let /C(-) be a kernel satisfying 


|1 — IC(x)\ < A|x| s , \/x £ JR \ {0} 


(28) 


with some A > 0. 
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(i) Assume that the density of A^ belongs to P(/3,L) with some [3 £ 1R and 
L > 0, and moreover 

a(j ) < e~j a , j e N, for some a* > 0. 

Then under the choice V n = 7 j 1 /( 2 ' 3 + 2s + 3 ) ; the MISE of the estimator v n 
is bounded as follows: 

MISE(y n ) < n - 2 «/(2/3+ 2 ^+3) ; „ oo. 


(ii) If the density of A a0 belongs to the class £(a,L) and 

a(j) < j~ a , j £ N, for some a* > 2, 
then under the choice 

1 s + 2 

V n = — log(n) - --- log(log(n)) 

la a 

we have 

MISE(D n ) < log _2s (n), n —> oo. 
Proof. Proof is given in Section 6.2. 


;□ 


The next theorem shows that the rates obtained in the previous theorem are 
optimal up to a logarithmic factor. 

Theorem 4.6 (lower bounds for v n ). Fix some seNU {0}, R > 0, a > 0, 
/3 > 0, L > 0 and define 


ip n {ir) := (p n (ir, p) = 


l s/(2/3+2s+3) log-P(n), if 7T € U(/3, L), 


log s (n), 


j/tt e £(a, L), 


for any p > 0 and any probability density tt £ "P(/3, L ) U£(a, L ), Then for some 
p* > 0, it holds 


inf sup 

TeS(s.-R) 

TT T eV(P,L)UE(a,L) 




<gm 

r 



— v(x)\ 2 dx 


> 0, (29) 


where the infimum is taken over all possible estimates v n of the function v based 
on i.i.d. sample Xi,... ,X n from the distribution ttj - of A^ := J 0 °° e - ^* dt such 
that the Levy triplet T of (ft) belongs to G(s,R). 

Proof. Proof is given in Section 6.3. D 

An important condition of Theorems 4.4 and 4.5 is (26), which means that 
the sequence Xq, Xi,..., X n is exponentially a-nrixing. Since /3-mixing coeffi¬ 
cient between two sigma-algebras is larger than or equal to the corresponding 
a-nrixing coefficient, it is sufficient to show that Xq. X-\..... X n is an expo¬ 
nentially /3-mixing sequence (see Section 1.1 from [8]). For the case of the GOU 
processes (1), the latter question was addressed in [12]. The sufficient conditions 
for exponential /3-mixing given in [12] are: 
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1. the distribution of A x has a Pareto-like asymptotic behaviour, that is, 

P{^4 00 > x} x Cx~ a as x — > oo 
with some a > 0 and C > 0; 

2. there exist A > 0, B > A and h > 0 such that ip(A) = 0 < oo with 
%l> given in (7), and 


E 


, — f.h 


e^ u ~du 


< oo. 


As it is proved in [22], both conditions are guaranteed by the positiveness of /i 
and the existence of a positive zero of the function ijj(-). We refer also to [21] 
for some further results in this direction. 


5. Simulation study 


Example 1. Consider the subordinator r t with the Levy density 

v(x) = abexp{—bx} E{x > 0}, a,b > 0. (30) 

Note that in this case, A = v{u)du = a. Define a Levy process 

= /A + crW t + Tf , (31) 

where W t is a Brownian motion. The Laplace exponent of is given by 

Hz) = * \ a2z + 5 ^) • ( 32 ^ 

In [9], it is shown that the exponential functional A = J 0 °° e - ^* dt is finite for 
any /x and a , and moreover the density function 7r of A 0 0 satisfies the following 
differential equation 


- —x tt"(x) + 


- (3 — b) + fi ) x — 1 


t'(x) 


(1 ~ b) ( + A 1 ) - a + ~ 


7 t(x) = 0. (33) 

\ 4 J X J 

Some special cases are considered below: 

1. In the case n = 0, a = 0 (pure jump process), this equation has a solution 

7T!(a;) = Cx b e~ ax I{x > 0}, (34) 


and therefore Aao = G(b + l,a), where G(a,/3) is a Gamma distribution 
with shape parameter a and rate j3. 










D.Belomestny and V.Panov/Statistical inference for GOU 


15 


Re 


Im 
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lm(s) 


lm(s) 




Fig 1. Plots of theoretical (blue dashed) and empirical (red solid) Laplace exponents in Ex¬ 
ample 1. Real, imaginary parts and absolute values are presented. 


2. If /i > 0, a = 0 (pure jump process with drift), then 

tt 2 (x) = Cx b ( 1 - !{0<x< 1 hi}. (35) 


In this situation = B(b+ 1 ,a/n)/n, where B(a,f3) is a Beta - distri¬ 
bution. 

3. In the case fi 0, a ^ 0, the equation (33) also allows for the closed 
form solution. Assuming for simplicity a 2 12 = 1, fi = —(6+ 1), we get the 
solution of (33) in the following form: 


„,(x) = C exp { ± J h Q-), 


(36) 


where we denote by 1^ the modified Bessel function of the first kind, n = 
yja + 1/4, and the constant C is later chosen to guarantee the condition 
/ 0 °° 7 r 3 (x)dx = 1. 


For our numerical study, we assume that the data are generated from the dis¬ 
tribution of (2), where the process (£ t ) is defined by (31) with fj, = 1.8, a = 0, 
and the subordinator r t in the form (30) with a = 0.7, b = 0.2. A sample from 
the distribution of the integral A a0 can be simulated from the corresponding 
Beta-distribution, see (35). In the first step, we estimate the Mellin transform 
A i(z) for z = u + in with u = u° = 29 and u = u° + 1 = 30 and v lying 
on the equidistant grid between —30 and 30. Next, we estimate the Laplace 
exponent of £ by the formula (12). Figure 1 graphically compares the pro¬ 
posed estimator of the Laplace exponent </>(u° + in) with its theoretical values 
(n + a/(b + u° + in)) • (u° + in) . 
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Fig 2. Boxplots for the estimates of ji = c and A = a for different sample sizes n based on 25 
simulation runs. 


Estimates for the parameters ft and A = a are given in (15) and (17), respec¬ 
tively. The boxplots of this estimates based on 25 simulation runs are presented 
on Figure 2. 

Example 2. Consider the compound Poisson process 


6 = - log q ri^J , 

where q £ (0,1) is fixed, N t is a Poisson process with intensity A and rjk are 
i.i.d. random variables with a distribution L. The integral A aa admits the rep¬ 
resentation 


Ax> = / q gt dt = y q Sn (T n+ 1 - T n ), 

J 0 n—0 

where T n is the jump time of N, i.e., T n = inf {t : N t = n}, and S n = Y^k= l Vk- 
Note that if % take only positive values, then is a subordinator. For the 
overview of the properties of the integral A aQ in the particular case rjk = 1 (that 
is, is a Poisson process up to a constant), we refer to [7]. 

Fix some positive a and consider the case when L is the standard normal 
distribution truncated on the interval (a, +oo). The density function of L is 
given by 

Pc(x) =p{x)/(l - F(a)), 

where p(-) and F(-) are the density and the distribution functions of the standard 
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Abs 



Fig 3. Plots of theoretical (blue dashed) and empirical (red solid) Laplace exponents for 
Example 2. Graphs present real, imaginary and absolute values. For v E [—3,3] the curves 
are visually indistinguishable. 


Normal distribution. In this case, the Laplace exponent of is equal to 

\2 ,2 ' 


= A 


l-F(a + (logq)z) J (logg) r 2 
1 - TT,qS)-“P-5- 


where the function F{-) can be calculated for complex arguments from the error 
function: 

F(z) := ^ ^erf + 1 ^ , where erf(^) = ~^= J e~ s ds. 


In this example, we aim to estimate the Levy measure of the process (£ t ), which 
is given by 


v(dx ) 


-- —— p{x)I{x > a}dx. 

1 — r (a) 


For our numerical study, we take q = 0.5, a = 0.1, and A = 1. First, we estimate 
the Laplace exponent by (12). The quality of the corresponding estimate at the 
complex points z = u° + iv with u° = 1 and v £ [—5, 5] can be visually seen 
in Figure 3. Next, we proceed with the estimation of the Fourier transform 
of the measure D(x) := e~ u X v(x ) by applying (18). For the last step of the 
Algorithm 2, i.e. the reconstruction of the Levy measure by (19), we follow [3] 
and use the so-called flat-top kernel, which is defined as follows: 


/C( x) 



b — 1/(1* I —0-05) 
1 - 1*1 


|cc| < 0.05, 

0.05 < |a:| < 1, 

M > 1. 


The quality of the resulting estimate v n is shown in Figure 4. 
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Fig 4. Left: plots of the Levy density (blue dashed line) and its estimate u n (z) (red solid). 
Right: the imaginary part of the estimate i'n(z) (red solid) and the line Y = 0 (blue dashed 
line). 


6. Proofs 

6.1 . Upper bounds for the quadratic risks of /i n and 

The next proposition is the main technical result for this section. 

Proposition 6.1. Let £* be a Levy triplet from Q(s,R). Suppose that the se¬ 
quence Xq, Xi, . .., X n of observations of the exponential functional Aqq := 
liniT->.oo At = f 0 e~& dt is a-mixing and strictly stationary. Denote the mixing 
coefficients of the sequence Xo, Xi,..., X n by o(s). 

Then for any p G {0,1,..., n} we have 
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provided, Y^jLi al £ (j) < 00 f or some e > 0 and the sequence V n satisfies 


sup 


= o(n 1 ^ 2 ). 


ve[o,v„] I M.{u° + iv + 1)| 

Proof. 1. Denote Y(z) := <j>(z) = z ■ M.(z)/M(z + 1), then 

POO /»oo 

M= / W/i,n(«)Im [h(«° + iv)] + / ■u; Ati „(w)Im[J r [p](-u)] du 

•to Jo 


(37) 


and we have 


nOO 

Hn- n = / W/i,n(v)Im [F n (v° + iv) - Y (u° + iv)] dv 

Jo 

nOO 

~ / w Mjr! (u)Im[J r [^](-t;)] dv 


= Im 


w^ n (v)S n (u° + iu) dv 


- Im[£>„(u°)] 


with 


S n (u° + iv) = Y n (u° + iv) -Y(u°+ iv), D n (u°) 


i{v)F[v](-v)dv. 


Note that 


(Pn - /^) 2 

< 2 • E 

Im 

■ /»oo 

/ w M , n (v)S n (v° + iv) dv 

)1 

+ 2\D n (u°)\ 2 

L J 


_V 

Jo 




2. Since 

S n (z) _ M„(z) _ M(z) 

z M n (z + 1) M(z + 1) 

M n (z)M(z + 1) -M(z)M n (z + 1) 

M n (z + 1)A 4(z + 1) 

[Mn{z) - M(z)\ Mn(z + 1 ) - [Mn(z + 1) - M(z + 1)] M n (z) 

+ l)M(z + 1) 

[Mn(z) - M(z)\ Y n (z) [Mnjz + 1 )-M(z + 1)] 

M(z + 1) 2 M(z + 1) 

[Mn(z) - M(z)] _ S n (z) [Mn(z + 1) - M(z + 1)] 

M(z + 1) 0 M(z + 1) 

Y(z) [M n (z + 1)-M{z + 1)] 
z M(z + 1) 

we get 

Sn ■ (1 + 7?2,n) = — ' 7?2,n + 
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with 

p M JM n (z) - M(z)] [M n {z + l)-M{z + lj\ 

Rl -" ( - ) = " M(z + 1) ’ fl2 -” W = -M^TT)-' 

Following the lines of the proof of Theorem 1.5 from [ 8 ], we get 

1 


E 


I M n {z) — M(z)Y 


n* 


E cov (x r \*n 


0<fc,j<n—1 


= ^ Var (X 5 - 1 ) + ^ E (l - l) Cov (A/T 1 ,^- 1 ) • (38) 

1 fc=i ' ' 

Note that the sum in the last representation converges as n —>• oo, because by 
Davydov’s inequality 


jCov(A 0 z - 1 ,Ar 1 )| < ^ (2 a(k)f- 2)/r (fi [x,f " 1)r ])'E (39) 
and therefore the series J2 (Aq~ : , A^ _1 ) is convergent if r = 2 /e. 


We have E 
result 


\A4 n (u° + iv) — A4(u° + iu)| < n 1 uniformly in v £ K. As 

1 


E 


\R2,n{u° +m)T 


< 


n • | M(u° + iv + 1)| 2 


The condition (37) implies now that sup^rgyp \R 2 , n {u° + iu)| = op(l). Fur¬ 
thermore, we have 


Var 


Rl,n(u° + i v)w tlt n(v) dv 

Cov(M n (u° + iui), Mn(u° + 1V 2 )) 


pOO poo 


/ 0 Jo M{u° + iv\ + l)M(u° + i ^2 + 1) 

• (u° + m)(u° - iv 2 )w^ n {v{) w^„(v 2 ) dvidv 2 - 


Similar to (38), we consider a representation 
Cov(M n {u° + mi), M n (u° + iv 2 )) 


1 


n —1 


go(vi,v 2 )+ 2^2g j (v 1 ,v 2 )+ 2 E 9j(vi,v 2 ) 


i =1 


where gj(v\,v 2 ) := (1 - j/n) • Cov [xf +1Vl l t xf +1V2 1 
Applying once more Davydov’s inequality, we get 


j=p +1 

, j = 0 ..(n - 1 ). 


( 2a C?)) (r 2)/r ( E 


A, 


(u°-l)r 


1 \ 2 / 


X z/r 

) ’ 


(40) 
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Now using the Cauchy-Schwarz inequality we get 


Var 


LJo 


Rl,n(u° + i V)w^ n (v) dv 


< 


\u° +iv\ 2 |w M ,„(e)r 


/ o \M(u° + iv + 1)|^ 


dv 


- 

j=p +1 


|W° + iv\ \w^ n (v)\ 


/ o \M(u° + iv + 1)1 

Finally using the fact sup^g^vy |i? 2 ,n(w° +iv)|^ = op(l), we derive 


dv 


Var 


S n {u° +w)w fl ,n(v) dv 


< 


r \u°+w\ 2 \w^ n (v)f 

la \M (u° + iv + 1)| 2 


dv 


+ 

j=p+1 


00 |u° + iv| |w jLt) „(«)| ^ 2 


Jo \M(u° + iv + 1)| 
3. Turn now to the term D n . By the Plancherel’s identity 


n(u) .. . . s r . 

[(—re) fF[v\(—v)\ dv 


(—iu) s 
Wp,n{v) 
(—iu) s 




dv 


= 2tt 


i/ ( s) (x) T~ x [w^n(-)/ (- i-) s ]{x)dx 


< 2 7 r||^)|| 00 ||J- 1 [ WMirj (.)/H.) s ]||L 1 . 


□ 


Proof of Theorem 4.4 


(i) Suppose that n £ V(f3,L) and a(j) < e ■ 7 “ , then by taking p = clog(rc) 
for c large enough, we arrive at 


E 


I Bn - 


< 


V n 4 log(n) 




M 2/3+2 kp(u/V'n )| 2 d« + TC 2(s+2) 


< n - 1 log(n) Vn 20 - 1 + V- 2(s+2) , 


E 


|A„ - A| 2 ] < Vn l0g(n) [ “ \v\^ \w x (v/V n )| 2 de + V r 

J n Jo 

< n- 1 log (n)V n 2f3+1 + V- 2 < s+1) 

By taking V n = n V(2/3+2s+3), we get 
E 


— 2(s+l) 


and 


E 


iMn - Ml 2 ] < n _2 ( s+2 )/( 2 ^+ 2s + 3 ) log(n) 

|A„ - A| 2 ] < n -2( S +D/(2/3+2 S +3) Jog(n) _ 
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(ii) Suppose that n £ £(a, L), then by taking p = 0, we get 


E 


E 


IMn - 


|A„-Ar 


< 


K 


-4 


f Vn \u° + in| 1 w,j,(v/V n )\ 
I o exp(—a|v|) 


dv 


i y-2(s+2) 

i v n. 


1 


< - exp(2aU„) + V n 2(s+2) , 


< 


n 
V~ 2 

* n 


n 

V 2 

v n 

n 


| U° + in| \w\(v/V n )\ 


dv 


Jo exp(-aH) 

< — exp( 2 al 4 .) + V~ 2 ( S+1 J. 

Tl 

Under the choice V n = ^ log(n) — log(log(n)), one derives 
E 


T/ — 2(s + l) 
1 v n. 


and 


E 


iMn-Ml 2 ^log 2(s+2 \n) 
|An-A| 2 l <log- 2 (s+ 1 ) (n). 


6.2. Upper bounds for MISE{u n ) 

Proposition 6.2. Let the assumptions of the Proposition 6.1 be fulfilled and let 
the kernel JC(-) satisfy the assumption (28). Then the mean integrated squared 
error of the estimator D n (x) satisfies the following asymptotic relation 


MISE(D n ) < - [ 
n Jm 


1 /' \u° +iv\ 2 [JC(v/V n )} 2 


ir | M(u° + iv + l)r 


dv 


+CiU„ 3 -E 


(Mn - hY 


+ C 2 v n ■ E 


(a„ - a y 


c, 


AL 

V 2s 


with some C\, C 2 , C 3 > 0. 
Proof. Recall that 
1 


v n {x) = / e ivx F[v}(- V )JC(-v/V n )dv = E- 1 [T P (-)IC(-/V n )](x), 

^ J IR 


and 


/[^](-u) = -Y n (u° + iv) + p n ■ (u° + iv) + \ n , 

F[v] (— v) = —Y(u° + i v) + p • ( u° + iv) + A. 
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By the Parsenval’s identity, 


MISE = -'-E 

2tt 


= 

2tt 


IM 


\T[v n ](v) - F[v\{v)\ dv 


i JR. 


F[v]{v)K(v/V n ) - F[v]{v) 


dv 


-E 


2tt 

< -E 

7T 


/JR 


(P[D]{v) - T[u]{v)) K{v/V n ) + (JC(v/V n ) - l)T[v]{v) 
(t[v]{v) - f[V]{ vj)K(v/V n ) 


dv 


+ -E 
7 r 


/JR. 


dv 


\(JC(v/V n )~l)F[D}(v)\ 2 dv 


3 . T x 1 

< — (Ji + Ji + J 3 ) H—<^4j 

7T 7T 


where 

Ji := E 
J‘2 ■= A n • E 
J 3 := B n ■ E 


[ | Y n (u° + iv) - Y{u° + iv )| 2 [K{v/V n )} 2 dv 
J JR 


(Mn - //)" 


with := / |u° + iz ;| 2 • [/C(v/V^)]“ dv, 
JlR 


(A„ - A ) 2 with := / [/C(v/14 )] 2 dv, 


/JR 


J 4 := 


/JR 


|(/C(v/E n ) - 1 ) J"[/j](v)| dv. 


The treatment of Ji is based on the observation that 


Y n (z) - Y(z) >c R 1>n = z 


[M n {z) - M{z)\ 
M(z + 1) ' 


We get that 


Ji x 



\M n (u° + iv) — M(u° + i-u) | 2 


|jj° + iv | 2 [/C( v/Vn)] 2 dv 
| M(u° + iv + 1)| 2 


As it was shown before, E 
Therefore, 


\M n {u° + iv) - M(u° + iv )| 2 


< n 1 , see (38)-(40). 


T < I . [ \u°+iv\ 2 {IC(v/V n )} 2 

~ n Jm |A4(w° + iv + 1)| 2 

To complete the proof, it is sufficient to note that 


A n ~V£-[ y 2 [K.(y)} 2 dy, 

JIR 


B n = V n ‘ [ [)C(y)} 2 dy, 
J ]R 
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and 


J 4 < A 



2s 


\T[v\{v)\ 2 dv < 


AL 

V^‘ 

v n 


□ 


Proof of Theorem 4.5 

(i) Recall that if n £ V(f3,L), then 


< n-H og{n)V n 2 ^- 1 +V- 2 ^ s+2 \ 
|A„-A| 2 < n- 1 ! og(n)V^ +1 +V- 2 ( a+1 >, 


see the proof of Theorem 4.4. Taking into account that J\ < n 1 V^ /3+3 , 
we arrive at 

MISE(i5 n ) < n~ 1 V 3+2 ^ + n" 1 log {n)V 20+1 + V~ 2(s+1) + V~ 2s . 
Choosing V n = n 1 /( 2 ^ +2s+3 \ we get 

n~ 1 V 3+2 P = V~ 2s > P- 2(s+1) , 

and therefore 

MISE(f'„) < n~ l log (n)V 2/3+1 + V~ 2s < n ~^/. 

(ii) Similarly, we derive the upper bound for the class £ (a, L ). Recall that 


E 

E 

and therefore 


\Bu-b\ 2 ^ n 1 exp(2 aV n ) + V n 2( ' s+2 \ 

I An A | 2 < n“ V„ 2 exp(2aP„) + F“ 2(s+1) 


MISE(P n ) < n- 1 log (n)V°e iaV " + V~ 2s < (logn) 


6.3. Lower bounds for MISE 

Proof of Theorem 4.6 . The general idea of the proof is to apply Theorem 2.7 
from [32]. This theorem yields that (29) holds, if there exists a parameterized 
set of Levy triplets 


T e = (l,0,^)c£(s,i?), d€{ 0,1} L 

for some s £ NU0,1? > 0, i > 0 and a set of parameters {0^\ j = 0,..., M} 
such that the following two properties hold. 
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(i) For any 0 < j < k < M, 

/ Wqu){x) - v g(k ){x)\ 2 dx>2ip n 

J R 


(41) 


(ii) Denote by ng.. j = 0,..., M, the probability distribution of the exponential 
Levy model Aj j00 = / 0 °° e~^ j ' s ds : where £y s is a Levy subordinator with 
triplet Tg r Then 


Yi X! K ^ eu) . ) < ^ io g( M ). 

3 =1 


M 


(42) 


for n large enough, where I\ stands for the Kullback-Leibler divergence 
between models, and x £ (0,1/8). 

Below we present a detailed proof for the polynomial case. 

1. Presentation of the models. Consider an exponential Levy model 
4lo,oo = / 0 °° e~^°' B ds, where £o,s is a Levy subordinator with a triplet (1,0, Pq) 
and Vq{x) = abe~ bx for some 0 < a < 1, 0 < b < 1. It is clear that (1,0, v 0 ) £ 
1/(0,-R) for some R > 0 and the Laplace exponent of £ 0) s is given by 


/•OO 

cf>o(z) =z+ (1 — e~ xz ) vq(x) dx = z 

Jo 


1 + 


Re(z) > -b, 


see Example 1 from Section 5. For the case of general classes Q(s, R) with s > 0, 
we could take a Levy density of the form vo(x) = b 1+s x s e~ bx /T(s + 1). 

Fix some L > 0 and let us construct now a parameterized set of Levy triplets 
Tg = (1,0, vg), 8 £ {0,1} L , with Levy measure ug defined by 

v e (x) \= v 0 (x) + 6 ■ A e (x), 

where S > 0 small enough, 

A g (x) := (g e (x) + a(gg-kexp(-b-)(x)))' : 

2 L 

g e (x) ■■= ^ 8k-LCOs(kj L x)g 0 {x ), 

k=L -\-1 

dk~L stands for the {k — L)-th component of the vector 8, —> oo as L —> oo, 

and 

go(x) := x~ 3 ^ 2 exp(— 1/x), x > 0. 

2. Distributional properties of the models. In this step, we perform 
some technical calculations, which will be used later. It holds 


noo 

C[/\ e \(z)= / e~ zx A e (x)dx 
Jo 


= z 


1 + 


z + b 


e zx gg(x)dx 


= Mz) -£[ge\(z), 
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where C[gg\(z) is the Laplace transform of the function go(-), which is equal to 
1 2i 

£[9e](z) = j E dk ~ L [ £ [^o K z + i iLk) + C[g Q \(z - b/ L k)\. 

k=L -\-1 

We see that / 0 °° A g ( x ) dx = 0 and 

<j>e(z) - </>o(z) = 6(j> 0 (z)C[gg](z), 

where is the Laplace exponent of a Levy process with the Levy triplet 
7 e- Furthermore, the Laplace transform of go is given by 

£[go\(u + iv) = s/i Te~ 2{z + +lz -'> 

with 2 z± = Vu 2 + v 2 iu. The Mellin transform of the density irg corresponding 
to the Levy model Tg satisfies the following functional equation 

Mg(z) _ 4>g{z) Mg(z + 1) 

M 0 (z) <j>o(z) Mo(z + 1)' 

Since 

W >~ 1=i£[wl(2) ’ 

and 

OO 

^2\£[gg\{z + k)\ < C exp (-^/2Re(z) - \J 2 |Im( 2 )|) 

k =1 

oo 2 L 

•E E exp(- a/ 27 L j - V2k) 

k=lj=L+l 

< C' exp ^/2Re(z) — sj2 |Im(^)|^ , Re(z) > 0, 
we derive the following infinite product representation for the ratio A 'lg(z) / A4q(z ) 


Mg(z) 

M 0 (z) 


tl(l + ££[3e](z + fc))- 

k =0 


Furthermore, it can be proved that 


Mg{u + w) 
Mo{u + in) 


< c5|£[g e ](u + in)| 


for some absolute constant c > 0. Note that the random variables Ag j00 = 
e~^ B ' s ds with being a Levy process with the triplet Tg , satisfies 0 < 
Ag,oo < 1 a.s. Moreover the density po of the r.v. Ao i00 has the form 


TT 0 (x) 


1 

B(b — 1, a) 


X b {l — X) a 1 l{o<a:<l} 
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and the Mellin transform A4 0 (z) of A 0iOO is given by 


M 0 (z) 


B(z + b, a) 
B{b — 1, a) 1 


(43) 


see Example 4.3. 

3. Class Q(s,R). In this step, we check that constructed models Tgu), j = 
1 ,M belong to class Q(s, R) with s = 0 and some R > 0. We have for any 
9£{0,1} l , 


( \v\ 2s \T[v e ](v)\ 2 dx < [ \v\ 2s \F[v 0 ]{v)\ 2 dv 

Jr Jr 

+ [ \v\ 2s \F[v e \ {v)-T[v 0 \ (v)\ 2 dv 
JR 

< f \v\ 2s \f[v 0 ](v)\ 2 dv + 6 2 f \v\ 2s \B[A e ](v)\ 2 dv. 

Jr Jr 


The inequality |</>o(— m)| < c • |u| for v £ K, where c = 1 + a/b, implies 

[ \v\ 2s \T[A g ](v)\ 2 dv < c f |n| 2(s+1) \£[gg\(-iv)\ 2 dv 
Jr Jr 

= - [ Ui| 2 ( s + 1 ) . 


2 L 


C 

< - 

- 2 


d k-L (C[g 0 ](-iv + i j L k) 

k=L -\-1 

+£[go](~iv - i7i^))| 2 dv 

2 L r 

V / \v\ 2l ' s+1) \C[g 0 ]{-iv + i^y L k)\ 2 dv 

A —L-f-l 

2 L r 

+ | / M 2(s+1) |£[sf 0 ](-m-i7ifc)| 2 dv 

l _r i i JR 


k=L-\- 1 1 


+Rl > 


where 


Rl = 2^/ \v\ 4 C[go]{-iv-ij'y L )£[g 0 \(-iv-ik'y L )dv 

k^j Jk 

+2V [ \v\ 4 C[g 0 ](-iv+ ij'y L ')C[go](-iv+ ik'y L )dv 
, /, Jr 




It holds 


2 L 


\Rl\ < CLj2(nL) 2{s+1) eM-J^d) 
j =i 

< C'L 2(s+1)+2 7L 2(s+1) exp(- v / 2^) 
= 0 (l 2 ( s+1 ) +1 ) , 
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provided 7 ^ = clog 2 (L) for large enough c > 0. Hence f R |f| 2s |J r [A e ](v )| 2 dv is 
bounded if < 5 2 7 2 (s+ 1 ) L 2s + 3 = 0 ( 1 ). 

4. Upper bound for the L 2 -distance between elements of {ug}. 

Fix two vectors 9,9' € {0,1} L . We have 

I \vg(x) - vg>(x ) 2 dx = f^d 2 [ |</>o(-i v)C[gg - ge'](-iv)\ 2 dv 

JR 27r JR 

1 2 L 

= J2 {°k-L ~ 9' k _ L ) 2 


k—L -\-1 


2 L 


\</)o(-w)£[g 0 ](-w + i^Lk)] 2 dv 


27r 


<$ 2 e 


k—L -\-1 


\<j>o(-iv)C\go](-iv - i 7 z,k)| dv 


+ l s2R - 


where 


R l < 2 V [ \<j) 0 (-iv)\ 2 C[g 0 \(-iv - \j^ L )C[g Q \{-iv - ikj L )dv 

+2V [ \4>o{-'w)\ 2 C[g 0 \{-\.v + \j^ L )C[g 0 }{-iv + \k'y L )dv. 
,, Jr 


Consider, for example, 


+ ++| 2 dv = [ \4> 0 (-i{v+ 'y L k))C[g 0 ](-iv)\ 2 dv 

Jr 


= / iv+7z^r 

Jr 

= ilk 2 


b — +; + j L k) 
a 


b — i(v + 7 L k) 


V 2 < dv 

g-2%/ 2 H dn + 0(jLk). 


So we have 

2L 


E (° k - L - 9' k _ L ) 2 f \cj) 0 (-w)C[go\(-w+ i'y L k)\ 2 dv 

r , , •/ K 


/c—L+l 


2L 


2L 


— C7l E (9k -l — 9' k _ L ) k 2 + 


2 7 9 . 2 

O 7 l 


/c—L+l 


fc—L+l 


E (e k . L -9' k _ L ) 2 k 2 J 

L 

>c'tiL 2 j2m^9' k ), 


k =1 
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as L —y oo and p(6 , O') = J2k =1 7 ^ #&) > 0* Analogously, 


2L « 

X! ( d k-L-0' k _ L ) 2 \M-iv)£[go\(.-iv-i'y L k)\ 2 dv = C"^ l L 2 p{6,0'). 

k=L +1 

Furtheremore, one shows (see above) that 

\Rl\ = o(L 3 ). 

5. Choice of 9^°\ ..., 0^ M \ 

Our choice is based on the well-known Varshamov-Gilbert bound (see [32], 
Lemma 2.9), which implies that there are M > 2 L / 8 vectors £ 

{0, \ } L such that 

p(6» (j ), 6>( fe) ) > L/8. 

6. Upper bound for K(no,ng). 

By Parseval identity for Mellin transforms, we get 


K{ n 0 ,Tr g ) = 


rl We(x) - 7r 0 (a;)| J 


i n 0 (x) 

= f x~ b (l — x) 1_a \TTg(x) — 7r 0 (a;)| 2 dx 

J o 

< j x~ b \ r Kg{x) —-Kq{x)\ 2 dx 

Jo 

= = / \Mg((l — b)/2 + iv) — Mq(( 1 - b)/2 + iu)| 2 dv 

J — OO 

o2 r oo 

< \M 0 ((l-b)/2 + iv)\ 2 \C[gg\(u + w)\ 2 dv. 

J — OO 


So we get 


C A 2 r 

K{n 0 ,-rrg) < — 55 / |7W 0 ((1 - 6)/2 + iv)| 2 |£[5 0 ](-iu + i'yz.A;)| 2 du, 

l _r i i J IK 


/c=L+l 
v2 2L 


c 2 ^ J-/ /» 

+ ^- 51 / |A^o((l-6)/2 + iu)| 2 |/:[5o](-i^-i7L*)| 2 ^ 

2?r r , JR 


fc=L+l 


where 


Rl < 2 55 f \Mq{( 1 - b)/2 + iv)\ 2 C[g 0 ](-iv - ij^ L )jC[g 0 }(-iv - ik"f L ) dv 

W Jr 

+ 2 55 f \- M o((l~b)/2 + w)\ 2 C[g 0 ](-w + ij'y L )£[go](-iv + ik'y L )dv. 
, .Jr 
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The equation (43) implies that Mq(z) is finite for all z with Re(z) > 0 and 

T(u + i v + b) 


Mo(u + iv) = C(a,b) 
x C(a,b} 


T(u + i v + b + a) 

— a log(u+ii;+b) 


= C(a,b)({u + b) 2 +v 2 ) a/2 e iAl ^ u+iv +V , u 2 +v 2 


Hence 


\Mo(u + in)| >c C(a, b) ((u + b) 2 + v 2 ) 


2 , 2\~ a / 2 2 i 2 . 

1 1 , U + V —» OO 


and the density 7r 0 of + j00 belongs to the class V(a,C(a,b)) (see also Exam¬ 
ple 4.3). We have 

2 l . 

V / \M 0 ({l-b)/2 + iv)\ 2 \jZ[g 0 \(-iv + i"f L k)\ 2 dv = 0(L~ 2a+1 ) 

_ T , J* 


and 


/c—L+l 


2L 


/c—L+l 

Hence 


E [ \M 0 {{l-b)/2 + w)\ 2 \C\g 0 \{-w-i'y L k)\ 2 dv = 0(L 2o+1 ). 
r , Jr 


J] E / ' r ( 7r °’ - n ^ L 2a lo s( M )’ o 


M 


(44) 


for large n. 

7. Choice of L. To complete the proof, we choose L such that the conditions 
(41) and (42) are fulfilled. First note that since our model belongs to the class 
Q(s,R), we can take jl = clog 2 (L) and S 2 = 7“ 2 ^ s+1 + _2s_3 • 0(1), see Step 3 
of the proof for details. Second, comparing (44) with (42), we fix x = nd 2 L~ 2a . 
This leads to the choice of L as the solution of the equation 

L 2a+2 S +3 log 4( s +l)(-£) = nQ ^ 

Combination of the results from Steps 4 and 5 yields the condition (41), because 

f \vg(x) - Ug>(x )\ 2 dx > OiJ+fL 3 

JR 

= C 2 (logL)- 4s L- 2s 
= C 3 (log L) 4S 2a +2°+ 3 n - 2 s/( 2 “+ 2s +3) 

for some 01 , 02,03 > 0 and L large enough. This observation completes the 
proof. 
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